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We compare the actual WMAP maps with artificial, purely statistical maps of the same harmonic 
content to argue that there are, with confidence level 99.7%, ring-type structures in the observed 
cosmic microwave background. 



1. Introduction 

Since its discovery the cosmic microwave background 
radiation (CMB) has been the subject of very intensive 
research. The CMB anisotropy provides invaluable infor- 
mation not only about the Universe in the epoch of last 
scattering but also about the earlier epochs. It follows 
from the lucky coincidence that before the last scatter- 
ing the Universe was radiation dominated, and in the 
presence of radiation pressure it could not evolve and all 
the structures essentially remained frozen. Therefore al- 
though the temperature at the last scattering was rather 
low (~ 1 eV), searching for structures of different types 
in the CMB anisotropy can possibly reveal facts about 
physics at energies far above those presently accessible 
in accelerators (^10 TeV). 

In this paper we compare the real CMB maps in 3 
frequency subbands Wl - W3 as measured by WMAP 
[l[ with artificially produced purely Gaussian maps of 
the same harmonic spectrum, looking for possible ring- 
type structures in the temperature distribution (see [2] 
for a theoretical motivation) . In this comparison we used 
the HEALPix code [3j to visualize and handle the maps 
but we used our own code to produce artificial maps. We 
find 2 statistically significant ring-type structures on the 
real maps (both happen to be on the southern Galactic 
hemisphere) that have no analogs on the artificial maps. 
We compare several characteristics of the real maps and 
the artificial maps. For all of these characteristics the 
real maps are very different from the artificially produced 
random maps. We applied a quantitative test to measure 
the difference and it allowed us to conclude that with 
confidence level exceeding 99.7% the ring- type structures 
we find on the CMB real maps are not the result of a 
statistical fluctuation. 

2. Creation of artificial maps 

Our algorithm for obtaining the maps with artificial 
CMB temperature distribution was as follows: 

1. we choose a maximal multipole number L (in the 
present paper we use L = 1000) 

2. for each multipole number I, < I < L, we choose 
a suitable number Ni, a positive integer, which is 
the number of spherical harmonics contributing to 
the temperature for this given I 



3. for all integers ki, 1 < ki < iVj, we assign a direction 
{Ok, , <i>ki ) = (arccos(7rar fcl ), 2iry kl ) where x kl and y kt 
are random numbers from the uniform distribution 
in the interval [0, 1] 

4. we take the coefficient C\ from the known real 
WMAP spectrum, and define a function 1} = 
Ti (9, <fi) on the sphere by: 



N, 



k, = l 



ft > iKcosMMAptf*,))) (!) 



where co{9, cj>, ( 
a given point 

{KM 



bk t ) is the spherical angle between 
■) on the sphere and the direction 



5. we repeat points 2)-4) for each I, < I < L 

The final formula for the temperature on the artificial 
map is: 
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It depends on our choice of the integers L and the iVj's. 
There is no preferred direction on the sky since for each 
I we choose randomly Ni directions {9k n (j)ki) (therefore 
this procedure has the advantage over averaging over m 
in the standard approaches). If we now perform the usual 
harmonic analysis of this map we get a set of Ci's that 
differ from the WMAP C; 's used in the procedure but the 
difference can be made as small as we wish with the choice 
of sufficiently large iVj's. The artificial maps produced in 
this way are purely statistical, Gaussian and, depending 
on the chosen Ni , can reproduce any prescribed spectrum 
of Ci with arbitrary accuracy. 

Since the spectrum of C; obtained from the real maps 
is (especially for small I) not precisely known we have 
produced artificial maps along the average, the upper 
and the lower C/ curve. For all of them we obtained 
essentially the same result. 

3. Looking for ring-type structures 

We have excluded from all the maps the Milky Way 
belt (±0.4 radians above and below the Galactic equa- 
tor) to avoid introducing structures caused either by the 
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Milky Way contribution or by the subtraction procedure 
usually applied to the CMB maps. 

The procedure for looking for the ring-type structure 
for any given real or artificial map consisted in the fol- 
lowing: 

1. a grid of points with HEALPix (k=6) parametriza- 
tion spreading over the entire sphere has been cre- 
ated 

2. a function 6 e ,7 fc (7): 

!-e(cos(7fc - e) - cos7 fc )~ 1 7fc - e < 7 < 7 fe 
e(cos(7fc) - cos(7fc + e)) _1 7fc < 7 < 7*; + 6 
otherwise 

was defined (this form is dictated by the area mea- 
sure on the sphere and the condition that a con- 
stant temperature should give no contribution to 
the integral below) 

3. for N (N w 7000) directions {6\ 0*) (i = l,...,N) 
belonging to the grid, for all widths of rings e 
(e = 0.02, 0.04 or 0.08 radians) and for all angu- 
lar "radii" j k = 0.07 + k ■ 0.01 (k = 1, . . . , 23) an 
integral 

/ (6 , Tfe) (0W) = J do'e e , 7fc ( w (^,0 l ,0'>'))7 1 (^^') 

was calculated (lo(0 1 , cf> 1 , 6', 4>') is the spherical an- 
gle between (#%</>*) and (9' ,<f)')). In this way for 
a given e we get a set of N ■ 23 « 160000 values 

4. for a given e we order these N ■ 23 values of 
J( e , 7fc )(# l , <p l ) (separately for positive and negative 
values), group them in a number n (n = 1500) 
of bins with centers Xi and for all artificial maps 
we form the normalized cumulative distribution 
Fm{ x i) ( m stands for the map number and i stands 
for the bin number) 

5. at each Xi we average over all F m {xi) to arrive at 
the average distribution F(xi) 

6. we use the procedure for comparison of cumulative 
distribution functions introduced in for all ar- 
tificial and 3 real maps Wl - W3 we calculate 

n 

^ = -^E^ ln ( 1 - Fa (^)) ( 2 ) 

i=l 

where a is a positive number (we used several num- 
bers, all very large, much larger than n) and di is 
the number of points in the ith bin. Then, sep- 
arately for positive and negative distributions, we 
check the hypothesis that a given map belongs to 
the distribution F(x) using the formula given in 
[H - the formula quoted there is suitable since it 
allows for quantitative comparison of distribution 
functions differing mostly at the tails (as is the case 
here). 



The shape of the function £)7fc ((9 l , </)*) defines what we 
mean by the ring-type structure - the maximal integral 
is when the temperature is of the same shape. Since we 
gather in the same set all values of 7*, the same direction 
can give several large values of the integral so these values 
can be correlated. The division into smaller number of 
bins n (n <C N ■ 23) is intended to partially compensate 
for these correlations. 

4. Results 

It is commonly taken for granted (with the notable ex- 
ception of @) that the Cosmic Microwave Background is 
purely statistical being produced by the quantum fluc- 
tuations usually assumed to have taken place during 
inflation (as the solution in De Sitter space suggests). 
Therefore it was very unexpected for us to find signif- 
icant differences (with confidence level 99.7%) between 
the WMAP results and artificial maps (with the same 
statistical properties as WMAP) that we have created. 
We believe that the procedure of creating artificial maps 
described above, being different from the one usually 
used, is interesting on its own. It gives purely statistical 
gaussian maps, and reproduces the prescribed harmonic 
content with the desired accuracy. 

We have produced many artificial maps (to generate 
an artificial map with our procedure is fast) but the pro- 
cedure for looking for ring-type structure required the 
number of integrals for any given map to be very large 
(N • 23 ?» 160000) and the CPU time needed was the 
main reason to limit ourselves in the present paper to 
100 artificial maps only. 

The differences between real and artificial maps were 
both qualitative and quantitative. First of all, for e = 
0.08 on frequency subband maps Wl - W3 we find two 
directions around which there is a significant concentra- 
tion of circular structures distinguished by the large val- 
ues of It e<JQ \{6,4>). The Galactic coordinates of these 

two directions are approximately (6i,<pi) = (2.6,3.7) 
(which correlates with the so cold "cold spot" p-0|) and 
(#2,^2) = (2.6,2.9) (this is new and has opposite values 
of the integral - the highest values from this direction 
are shown on Fig. 2). We interpret these two direc- 
tions as two centers of circular rings in the temperature 
distribution. The reason for this interpretation is as fol- 
lows: in a close vicinity (about 8 • 10~ 3 sr) of the first 
of these directions, there are about 80 centers of circular 
structures (and there are none outside). A similar situa- 
tion occurs around the second direction - there we find 
about 40 centers in about the same solid angle as in the 
first case. We also found many more potential centers of 
ring-type structures but since they were less significant 
statistically we do not include them in the analysis. It 
is perhaps worthwhile to note that the positions of both 
of the main two ring-type structures and the less signif- 
icant ones correlate approximately with the positions of 
the structures found in [8]. 

The quantitative comparison between the real Wl - 
W3 maps with the artificial maps regarding the existence 
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of structures with e = 0.08 yields the following: 

For the procedure described above we get for the real 
maps Wl - W3 with n = 1500, a/n = 2000 the values 
A = 107.9, 133.5, 84.8 respectively, while only 6 out of 
100 artificial maps have A > 1 (the largest value among 
these 6 maps is A ~ 52, the other 5 have much smaller 
A). According to [H the probability of A bigger than a 
for large a/n is approximately equal to (1 <C a <C a/n) 

, . s n , a 
P(A > a) = - In — (3 
a na 

so we get for the real map W2 the probability of a statis- 
tical fluctuation « 0.3%, which means that with 99.7% 
confidence level we can reject the hypothesis that the val- 
ues of the integrals on the real maps are purely statistical 
and state that the ring type structures on the background 
radiation maps are not statistical fluctuations. In favor 
of this statement we can also observe that for the value of 
/( £i7o )(#, (f)) greater than 3 on the map of frequency band 
W2 only 81 points are in the northern hemisphere and 
1430 in the southern one - the difference is similar for two 
other spectral types but there are no such differences for 



the artificial maps. 

It is interesting to note that for e = 0.02 and e = 0.04 
the differences between the real and artificial maps are 
much smaller. 

For comparison we have also imposed a condition that 
around the same direction (6, </>) we have both large posi- 
tive and large negative integrals. It turned out that none 
of the real maps shows any statistically significant pres- 
ence of such structures - it points to the presence of one 
sharp circular edge in the temperature distribution and 
not two or more (which may be there but milder). 

In figure 1 we have plotted the directions of the sus- 
pected ring-type structures on the real WMAP maps in 
3 frequency bands Wl - W3.(with the galactic disk ex- 
cluded, as described in the text). In fig 2. we have plotted 
the upper part of the histogram of the cumulative dis- 
tribution functions for all artificial maps and three real 
maps Wl - W3 (for positive value of the integral). 
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Fig.l. Position of centers of circular structures on the real maps Wl - W3 
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Fig. 2. The tail of the distributions of the integrals of the 100 artificial (blue) and 3 real (red) maps 



